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The problem of state selection when multiple metastable 
' states compete for occupation is considered for systems that 
are accelerated far from equilibrium. The dynamics of the 
supercurrent in a narrow superconducting ring under the in- 
fluence of an external electric field is used to illustrate the 
general phenomenology. 
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Many systems when driven far from equilibrium en- 
counter instabilities that lead towards new states or 
phases. Frequently there exist multiple states that can be 
selected following the onset of the instability. The deter- 
mination of the particular state that is selected is a com- 
plex problem of fundamental interest in a wide variety of 
fields . In addition to the complexity associated with 
the presence of multiple competing states, if the system 
is accelerated so that the environment evolves in time, 
then the state that is selected can depend in an impor- 
tant way on the driving force. The focus of this paper is 
on state-selection in "accelerated" systems when multiple 
metastable states compete for occupation. 

For the purposes of this paper an accelerated system 
is defined as one for which a control parameter is var- 
ied in time so that the system gradually progresses from 
stable to unstable regimes. For example, in a narrow 
superconducting ring under the influence of a con- 
stant electromotive force, the superconducting electrons 
are accelerated by the electric field and the supercurrent 
increases with time until the critical current is reached 
and the system becomes (Eckhaus) unstable. Similar be- 
havior could occur in direction solidification if the 
solidification cell is accelerated slowly, rather than pulled 
at a constant velocity, through a temperature gradient 
until the (MuUins-Sekerka) instability is encountered and 
the liquid/solid interface becomes unstable. In each of 
these scenarios the systems become unstable with respect 
to fiuctuations of certain wavelengths that lie within a 
band. When the size of the system is comparable to the 
length scales associated with the wavevectors in the un- 
stable band, the system is described as mesoscopic and 
the number of accessible states is finite. As illustrated 
in an extensive review by Cross and Hohenberg in- 
stabilities that result in this type of mesoscopic behavior 
are extremely common, occurring in many diverse fields, 
such as fluid dynamics, chemical reactions, material sci- 
ence and biology. 



The focus of this paper is on state selection in meso- 
scopic accelerated systems immediately following the on- 
set of the instability. The combination of the driving 
force that accelerates the system, and the mesoscopic 
system size that allows for multiple, isolated metastable 
states, leads to novel and interesting selection rules. It 
will be shown that the rate at which the system is driven 
through the instability plays a prominent role in deter- 
mining the probability that a particular metastable state 
is selected. As the decay is from states of marginal sta- 
bility, the selection is also influenced by the noise in the 
system. The dependences on both acceleration rate and 
noise strength are considered. 

For illustrative purposes consider the situation in 
which an infinitely long solenoid, carrying a current that 
increases linearly with time, passes through the center 
of a narrow superconducting ring of cross-sectional area 
A and circumference L = ^{T)i, where £,{T) is the tem- 
perature dependent correlation length. By Faraday's law 
of induction, a constant electromotive force (emf) V is 
induced in the superconductor, thereby accelerating the 
superconducting electrons. The dynamics of the (dimen- 
sionless) superconducting order parameter ip{x,t), where 
X is the longitudinal spatial coordinate and t is time, 
is described by the stochastic time-dependent Ginzburg- 
Landau equation 



(1) 



where to = TcLi'^sV/h) is a dimensionless measure of 
the strength of the induced emf. Throughout this paper 
the regime where w <C 1 is considered. In Eq. (|l|) tgl, 
is the Ginzburg-Landau relaxation time, and "0 satisfies 
the twisted-periodic boundary condition 'ijj(£ + x,t) = 
exp{iLL)t)'ip{x,t). This equation ignores the influence of 
the normal current and is valid for a voltage driven sys- 
tem in the limit of low normal state resistivity. The vari- 
able ry is a Gaussian random variable, with expectation 
values {7]{x,t)) = 0, and {r]{x,t)Tj*{x' ,t')) = 2D5{x - 
x')5{t — t'), where D is determined by the fluctuation- 
dissipation theorem |9|. 

For <C 1, the relevant current-carrying states of the 
superconductor are uniformly twisted plane wave solu- 
tions given by "0 = a/1 — exp{iqx), where q = mK -\- 
ijjtjl, and K = 2^/1. The dimensionless current density 
J of these states is given by J = {'ip*dx'ip — 'ijjdx'4'*)/'2i = 
q{l ~ q^). Thus the effect of the induced emf (which in- 
creases q linearly with time) is to wind the order param- 
eter, or equivalently, to accelerate the superconducting 
electrons. However, this acceleration cannot continue in- 
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definitely because J is a nonmonotonic function of q, and 
hence time, achieving a maximum value of Jc = 2/^/27 
a.t q ^ Qc ~ l/\/3- This saturation of the current at 
the critical current Jc, coincides with the loss of stabil- 
ity of states tj) dX q = qc- In other words, for q > qc, 
(PF{q)/dq'^ < 0, where F{q) ^ Fgl [■(/'] is the Ginzburg- 
Landau free energy of states -0- 

To understand the Eckhaus instability for finite size 
systems it is necessary to perform a linear stability analy- 
sis about the state ip, as the previous analysis only applies 
when £ = oo. Standard linear stability analysis gives one 
potentially positive eigenvalue that takes the form 

Xn{q) = -l + q^- kl + v/(l-92)2+Vfc2 . (2) 

The eigenvector associated with this eigenvalue is a linear 
combination of Fourier modes with wavevector q±kn and 
amplitude An, where fc„ = nK. The interesting feature 
of this eigenvalue is that it can become positive when 
q> Kl > qc, where = -^[1 + Thus, for 

finite size systems the instability is pushed to wavevectors 
greater than qc by an amount that depends on i In 
particular, for Km > q > ki, X„ is positive for all values 
of kn < nriK [|ll|. The dependence of A on fc„ is shown 
in the inset of Fig. |^ for several values of q. 

The growth of a single Fourier mode (with amplitude 
An) of wavevector q — kn, and simultaneous decay of 
Aq, corresponds to a decrease of the winding number 
W = (27r)~^ d(j){x)/dx, where (j) is the phase of -0, by 
an amount n. This phenomenon is known as a "phase- 
slip" as the total phase of the order parameter changes 
by an integral multiple of 27r. Physically, the supercur- 
rent decreases by a discrete amount when a phase-slip 
occurs. Phase-slip processes can also occur via thermal 
activation over an energy barrier and this process has re- 
ceived significant attention over the years . For 
D = 10""^, as long as w ^ 10"^'', the probability of a 
thermally activated phase slip occurring is exceedingly 
small 1^. Thus, unless the temperature is very close to 
the superconducting transition temperature Tc, where D 
is large, the system will almost always be driven to the 
Eckhaus instability before a thermally activated phase- 
slip can occur. Therefore, the transitions that are of 
concern in this work involve the decay from an unstable 
state, in contrast to previous work where the focus 
was on the decay from a metastable state. 

When ti> > 0, the system is driven to the point of insta- 
bility as the eigenvalues of each Fourier mode eventually 
become positive. As illustrated in Fig. ^ the n — 1 
mode becomes unstable first, then the n = 2 mode be- 
comes unstable, and so on. This implies that the system 
first becomes unstable with respect to single phase-slip 
processes, then double phase-slip processes, etc. If to is 
large enough, the n = 1 mode might not have time to 
grow to dominance by the time the n = 2 mode becomes 
unstable. This suggests that for small uj, single phase- 
slip processes should dominate the dynamics, but as m 
is increased there is a crossover to a regime in which 



double phase-slip processes dominate. As lo is increased 
further, double phase-slip processes should give way to 
triple phase-slip processes, and so on. 

The generic features displayed in Fig. ^ are common to 
many systems and come under the general classification 
scheme of Cross and Hohenberg as type Ilg. Thus, the 
dynamic competition between unstable modes discussed 
above is a phenomena that has relevance to many sys- 
tems. The precise determination of which of the modes 
will initially be selected following the onset of the insta- 
bility is a complex question that depends on the details 
of the individual system. For concreteness, the general 
problem of the dependence of the selected state as a func- 
tion of the driving force will be addressed for the super- 
conductor. In particular, the relative probability that a 
given state is selected will be determined by extensive 
computer simulations. Second, it will be shown that the 
qualitative features of these results can be understood by 
an analysis that is based on the properties of the growth 
rates A. 

To evaluate the probability of the occurrence of a given 
phase slip as a function of w, Eq. was numerically in- 
tegrated in time for a noise strength oi D ~ 10^'^ and 
a length corresponding to ni = £qc/{2i:) — 5 In 
Fig. eK, the probability of a type-n phase slip is plotted 
as a function of tj. As expected, for small lo, single phase- 
slips dominate. As lo increases further there is a crossover 
to a regime in which double phase-slips dominate. Fur- 
ther increase of lo results in a subsequent crossover to a 
regime in which triple phase-slips dominate, and so on. 

An example of the dynamics that lead to such results is 
shown in Fig. |[ In this figure, the winding number and 
current are plotted as a function of time for lo = bx 10"'*. 
This value of lo is in the crossover region between the 
single and double phase-slip dominated regimes. Clearly 
evident in this figure are the single and double phase slips 
in which W changes by one or two, respectively. Also seen 
in Fig. ^ are the discrete jumps of the supercurrent. 

As described earlier, the essential features shown in 
Fig. ||a can be understood using the properties of the 
growth rates A„. This idea can be made more concrete 
in the following way. Ignoring the nonlinear interactions 
between the different modes, the expectation of |^nP is 
given by ^ 

(|A„(t)|2) = ^e2-Wj[|*dt'e-2-(*'), (3) 

where cr„(t) = dt' Xn {q{t')), and angular brackets 
denote a noise average. After the onset of the in- 
stability the system evolves towards the fixed points 
^n = An exp{i{q — nK)x), where A„ = -^1 — — nK)'^. 
Eq. describes the initial evolution of the system after 
the Eckhaus boundary has been reached. In this non- 
interacting picture each amplitude (measured in units 
of An) can be thought of as an orthogonal coordinate 
in an n^-dimensional space. Thus, the natural measure 
of distance from the origin {An = 0) in this space is 



YZLi{\Anit)\'^) I Al- After onset of the Eckhaus instabil- 
ity, this sum increases rapidly and reaches unity at a time 
t* . Assuming that dXt = t* a phase-slip has occured with 
probability one it is natural to interpret (|A„(t*)p)/^^ 
as the relative probability of the occurrence of a type-n 
phase-slip. The probabilities calculated using this proce- 
dure are shown in Fig. Bp. 

It is clear from Fig. [Jthat the preceding analysis pro- 
vides a qualitatively accurate description of the state- 
selection probabilities, and their dependence on the driv- 
ing force Lo. Most notably, the values of lo at the peak 
positions agree very well with the numerical results. Nev- 
ertheless it is important to point out that the preceding 
analysis is only a plausible argument and is not system- 
atic. A quantitative description must include the subtle 
non-linear interactions that are an important element in 
determining state selection. Even at the present level of 
ignorance, however, the analysis presented here provides 
a qualitatively useful description of the state-selection 
probabilities, and their dependence on the driving force. 

The growth rates A are an extremely important fac- 
tor in determining the state selection probabilities. The 
preceding analysis accounts for these growth rates and 
therefore provides a qualitatively accurate description. 
The analysis also provides predictions for the dependence 
of P„ on the noise strength D, which may be more con- 
venient to vary in some experiments. Plotted in Fig. ^ 
are the probabilities of a type-rt phase slip as a function 
of -D, for a fixed value of uj, obtained from a numeri- 
cal simulation of Eq. (|l|). In Fig. the corresponding 
P„'s obtained from the growth-rate analysis are plotted 
for comparison. Once again, it is seen that the simple 
analysis provides an accurate qualitative picture. For 
the smallest values of D considered triple phase-slip pro- 
cesses dominate. This is because the time required for a 
given mode to grow to saturation diverges logarithmically 
as Z? ^ 0. Consequently, if D is very small, the mode 
amplitudes Ai and A2, for example, may still be very 
small by the time the growth rate of ^3 is significantly 
larger than the growth rates for Ai or ^2- 

One of the most interesting aspects of the phenomena 
exposed here is that the selection rules depend on both 
the intrinsic properties of the system and the external pa- 
rameters. To understand this connection more deeply, it 
is instructive to consider the characteristic growth times 
for individual modes. Typically, the characteristic time 
associated with the initial growth of an unstable mode 
is taken to be the inverse of the growth rate. How- 
ever, for accelerated systems the growth rate A starts 
out negative and passes through zero. Thus, |A^^| is not 
a relevant quantity as it diverges at the instability. To 
determine the characteristic time, consider Eq. (g|) for 
The quantity (T{t) achieves a local minimum 
&t t = tn = innjio so that a second order expansion 

about t„ yields (T„(i) w (T„(t„) -I- \ ^^^{t — t„)^, where 
AJj = <9A„/9(7|g=K^ . Inserting this expansion into Eq. (H) 
and assuming that w <C 1 gives 



(|A„(i)|2) =2Dr„£-iexp(z2(t))[erf(z„(i)) + l], (4) 

where z„(t) = (t — tn)/Tn and t„ — \JTJX^. The quan- 
tity r„ is the characteristic time for the growth of mode- 
n, and is interesting because it depends on the geometric 
mean of AJ^ and lo. Thus, the time scale t„ embodies 
in a natural way the importance of the combination of 
the intrinsic dynamics (A^) and the external driving force 
(c.). 

In summary, the problem of state selection in accel- 
erated systems has been shown to contain unique and 
rich phenomenology. Although the focus of this work 
has been on the dynamics of quasi-one-dimensional su- 
perconducting rings, the essential features of this partic- 
ular system are generic and should be observable in a 
diverse array of experimentally realizable situations. De- 
spite the success of the linear analysis, it is clear that new 
methods must be developed to explore this complex and 
important area of research in nonequilibrium statistical 
mechanics. Recent work on state selection in non- 
accelerated marginally stable systems suggests a possible 
systematic framework that could be extended to address 
the phenomena considered here. 
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FIG. 1. A as a function oi q = ut/e for k = K,2K,3K. 
Inset: A as a function of k for two values of g > ki such that 
the upper curve corresponds to the larger value of q. 

FIG. 2. State selection probabilities as a function of the 
driving force uj. Open squares, solid squares, open circles, 
solid circles and open triangles correspond to the probabil- 
ities Pi, P2, Pa, P4 and P5, respectively. Results of the nu- 
merical integration of Eq. (Q) and those of the linear analysis 
described in the text are shown in Figs, (a) and (b), respec- 
tively. 

FIG. 3. Dynamics of winding number (a) and supercurrent 
(b), for = 5 X 10"" and D = 10"^ 

FIG. 4. State selection probabilities as a function of the 
noise strength D. The symbols in this figure are identical 
to those in Fig. |^. Results of the numerical integration of 
Eq. (0) and those of the linear analysis described in the text 
are shown in Figs, (a) and (b), respectively. 



